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Abstract. We present deep, wide-field European VLBI 
Network (EVN) 1.6 GHz observations of the Hubble Deep 
Field (HDF) region with a resolution of 0.025 arcseconds. 
Above the 210 /iJy/beam (5a) detection level, the EVN 
clearly detects two radio sources in a field that encom- 
passes the HDF and part of the Hubble Flanking Fields 
(HFF). The sources detected are: VLA J123644+621133 
(a z=1.013, low-luminosity FR-I radio source located 
within the HDF itself) and VLA J123642+621331 (a dust 
enshrouded, optically faint, z=4.424 starburst system). A 
third radio source, J123646+621404, is detected at the 4a 
level. The VLBI detections of all three sources suggest that 
most of the radio emission of these particular sources (in- 
cluding the dusty starburst) is generated by an embedded 
AGN. 
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1. Introduction 

HST obser vations of the Hubble Deep Field North (here- 
after HDF, Williams et al. 1996 ) have provided deep, high 
resolution optical images of the distant universe. These 
have been vigorously followed up with complementary ra- 
dio, sub-mm, far infrared, near infrared and X-ray obser- 
vations (Ferguson et al. 2000 and references therein). The 
high sensitivity radio observations of Richards et al. 1998 
(hereafter R98), |Muxlow et al. 1999| (hereafter M99) 



Richards 2000 (hereafter R00) and Garrett et al. 2000 



have thrown new light on the nature of the faint microJy 
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radio source population. These observations, spanning a 
range of angular resolutions (0.15" to 15"), suggest that 
~ 60% of faint sub-mJy and microJy radio sources have 
radio luminosities, Li.4GHz < 10 23 W/Hz, steep spectra 
and angular sizes comparable or smaller than their parent 
galaxies. In the optical these systems are identified with 
moderate redshift, star- forming galaxies, often presenting 
peculiar morphologies and evidence for interactions with 
other galaxies in the field. The radio emission in these 
systems arises from a combination of young supernovae, 
supernova remnants (SNR) and, in particular, relic SNR 
emission. Approximately 20% of the faint radio source 
population arc identified with relatively low-luminosity 
AGN, with the remaining 20% identified with optically 
faint sources with no detections down to / = 25. 5 m in the 
HFF and I = 28.5™ in the HDF (R98, M99, R00). The 
optically faint systems are thought to be distant galax- 
ies, obscured by dust ( Richards et al. 1999| ). Such systems 
may only be detectable in the radio and sub-mm where 
the effects of extinction due to dust are less significant. 
Since the nature of the optical morphology and optical 
spectra of these sources is unknown, it is as yet unclear 
whether individual systems can be identified as pure star- 
burst galaxies, dusty AGN or perhaps some mixture of 
both phenomena. In principle, high resolution VLBI ob- 
servations can discriminate between these possibilities. 

In this paper, we present high resolution, VLBI obser- 
vations of faint radio sources in the HDF. The aim of these 
pilot observations, made with the European VLBI Net- 
work (EVN), was to probe for the first time, the milliarc- 
second (mas) radio structures of the faint micro- Jansky 
radio source population. Compared with the typical VLBI 
field of view, the HDF spans a relatively large area of sky 
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and was specifically chosen to be devoid of all strong radio 
sources (Si.4gh z < 2 mJy). The success of the observations 
described here relied heavily on the recently enhanced ca- 
pabilities of the EVN (in particular the introduction of the 
high sensitivity MklV data acquisition system), together 
with the use of the now standard techniques of phase- 



referencing (e.g. Beasley & Conway 1995) and wide- field 
VLBI imaging ( |Garrett et al. 1999| ). 



2. EVN Observations and Data Analysis 

Observations of the HDF region were made with the Eu- 
ropean VLBI Network on 12-14 November 1999 at 1.6 
GHz. A total of 32 hours observing time was split into 
two, 16 hour runs. The EVN included the 100-m Effels- 
berg (DE), 76-m Lovell (UK), 32-m Medicina (IT), 32- 
m Noto (IT), 32-m Cambridge (UK), and 25-m Onsala 
(SE) telescopes; recording data at a total sustained bit 
rate of 256 Mbits/sec (including both right (R) and left 
(L) hand circular polarisation). The Deep Space Network 
(DSN/NASA) 70-m telescope at Robledo (ES) was also 
used for approximately 11 hours of the total 32 hour run, 
providing data in one polarisation (L). 

The observations were made in phase-reference mode 
with a typical cycle time of 2.5 minutes on the phase- 
calibrator and 4.5 minutes on the HDF. The total on- 
source integration time on the HDF was ~ 14 hours. 
Two phase-calibrators were used, the primary calibrator 
J1241+602 (a compact S*i.4ghz ~ 455 mJy source ly- 
ing 2° from the HDF) and a fainter secondary calibrator 
J1234+615 (5i.4GHz ~ 30 mJy) lying only 23' from the 
HDF). This nearby, secondary calibrator was first detected 
in a targeted MERLIN survey of potential mJy calibrator 
sources located within 0.5° of the HDF. Having a nearby 
calibrator allowed us to reduce possible phase errors as- 
sociated with ionospheric disturbances (which may have 
been significant during this period of high solar activity). 
However, at the time the observations were conducted, the 
compactness of the secondary calibrator was not known on 
VLBI scales, it was thus considered prudent to only em- 
ploy it in every other tape pass, i.e. ~ 1/2 the total tape 
passes scheduled. 

The data were processed at the NRAO VLBA correla- 
tor in Socorro, NM, USA. The phase centre (and telescope 
pointing position) of the single pass correlation was coinci- 
dent with a 470 ^Jy radio source, VLA J123642+621331 
(R00), lying just outside the central HDF (in an adja- 
cent HFF). All subsequent data processing was performed 
with the NRAO MPS package. The visibility amplitudes 
were calibrated using the system temperatures and gain 
information provided by each telescope. The residual de- 
lays, fringe rates and antenna gains were determined from 
J1241+602, then interpolated and applied to the data. 
The secondary calibrator was found to have compact 
structure on EVN angular scales, and thus second order, 



phase-only corrections were also determined for those tape 
passes that included the secondary calibrator. 

We selected as possible targets, all sources located 
within 3.5' of the EVN phase-centre with MERLIN peak 
flux densities, Si.4GHz, m excess of 60 /zJy (as measured 
from uniformly weighted maps). Five target sources satis- 
fied these conditions, all of which turned out to be located 
within 2' of the phase-centre. A sixth, "blank" field was 
also chosen, in order to act as a check on the robustness 
of the imaging process. With the final gain corrections 
applied, the HDF data were Fourier transformed to the 
image plane and dirty images were simultaneously gener- 
ated for the six target fields. 

In order to image the entire HDF (an area of ~ 
5 arcmin 2 ), the HDF data were stored in the form in 
which the correlator generated them: 512 contiguous 0.125 
MHz spectral channels with a visibility averaging time of 
1 second. We estimate that for these data, the effects of 
bandwidth and time smearing are minimal, resulting in 
no more that a 10% reduction in source intensity for fields 
located more than 3' from the phase centre. A more im- 
portant limitation on the field-of-view is the fall off of 
the primary beam response of the larger telescopes in the 
array (primarily Effelsberg, Lovell and Robledo). The con- 
tribution made by these telescopes dominate the naturally 
weighted images presented here. The most distant target 
field presented in these observations is located 2 arcmin- 
utes from the telescope pointing centre. We estimate that 
our response is down by less than 15% in fields located 2 
arcminutes from the phase and telescope pointing centre. 

The naturally weighted dirty images are shown in 
Fig 1. The images were restored with a Gaussian ellip- 
tical beam (26 x 20 mas in PA= —31°). The rms noise 
levels range from 42 /iJy/beam (in the field with the 
strongest detection) to 33 /xJy/beam (in a field with no de- 
tection). This is higher than the predicted thermal noise 
(~ 23 /iJy/beam) which suggests we are currently lim- 
ited by errors in the calibrator phase solutions. With the 
maximum side-lobe level at 20% of the main dirty beam 
response, CLEANing the 2 fields with clear detections pro- 
duced no substantial improvement in the image quality. 

3. EVN Detections in the central HDF region 

At the 5<T detection level, the EVN observations clearly 
detect VLA J123642+621331 (at the phase centre) and 
VLA J123644+621133 (located within the HDF itself). 
Both sources are detected independently in both left and 
right-hand circular polarisation. The sources are separated 
by ~ 2'. There is a 4ct detection of J123646+621404, lo- 
cated within 20 mas of the measured MERLIN position. 
The EVN source positions, flux densities and sizes of all 
3 sources is presented in Table 1. The 1.4 and 8.4 GHz 
VLA flux densities (R98, R00) of the sources are also 
given. Two other target sources, J123646+621448 and 
J123652+621444, with peak MERLIN flux densities of 60 
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Fig. 1. EVN 1.4 GHz dirty images of the six MERLIN targets. The MERLIN source positions are identified by crosses. 
Contours are drawn at -3, 3, 4, 5, 6, 7, 8, 9 times the rms noise level ~ 40/iJy per beam. 

Table 1. EVN 1.4 GHz Detections in the HDF Region 



Source 


RA (±12 h 36 m ) 


DEC (+62°) 






Size 


Si. 4 


S8.4 




(<0 


(' ") 


/ijy 




mas 


/ijy 


/ijy 


VLA J123644+621133 (c) 


44.38626 ± 0.0001 


11 33.1767 ±0.001 


351 ±41 


351 ±41 


< 20 


1290 


477 


VLA J123644+621133 (j) 


44.38933 ± 0.001 


11 33.1162 ±0.002 


190 ± 41 


214 ± 77 


< 30 






VLA J123642+621331 


42.09078 ± 0.0002 


13 31.4254 ±0.001 


248 ± 37 


248 ± 37 


< 20 


467 


73 


VLA J123646+621404 


46.33177 ± 0.0002 


14 04.6893 ± 0.002 


180 ± 38 


180 ± 38 


< 20 


179 


168 



and 103 /xJy were not detected (> 3cr level, see Fig. 1). 
Note that since the primary calibrator (J1241±602) has 
an ICRF position, the absolute positions of the sources 
presented in table 1 are estimated to be accurate to bet- 
ter than 10 milliarcsecond. These are the most accurate 
ICRF positions associated with the HDF region. 



3.1. VLA J1236U+621133 

VLA J123644+621133 is identified with a bright (I ~ 
21.05" 1 ) z=1.013, red elliptical galaxy (R98). The radio 
source is clearly identified by its large scale radio struc- 
ture and luminosity as an FR-I (R98, M99). The EVN 
(see Fig. 1) detects a compact core component (c) with 
an estimated size < 20 mas, implying a brightness tem- 
perature, Tfc > 4 x 10 5 K. We also detect an additional 5a 
component (j) located ~ 60 mas south of the core compo- 
nent. Presumably this source (j) is a knot in the extended 



VLA/MERLIN jet that also trail southwards of the core in 
a similar position angle. Note that the MERLIN position 
for this source (see Fig. 1) is strongly affected by the ex- 
tended jet structure, lying mid-way between the compact 
core (c) and the knot in the jet (j). The VLBI structure 
observed in this source, supports the view (M99) that the 
rare "classical" radio sources detected in the HDF and 
HFF, are the mJy tail of the AGN radio source popu- 
lation. It is worth noting that VLA J123644±621133 is 
probably an order of magnitude more distant than any 
other classical FR-I radio source so far studied with VLBI. 



3.2. VLA J 123642+621331 

VLA J123642+621331 lies just outside the HDF, in an 
adjacent HFF. At radio wavelengths it is a relatively 
bright, steep spectrum radio source (Si.4GHz ~ 470 /iJy) 
with no optical counterpart to I < 25 m (R98). Deep 
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HST NICMOS 1.6 /im imaging clearly detects an ex- 
tremely red 23. 9 m counterpart to the radio source with 
a disk profile, and spectra obtained by the Keck II tele- 
scope show a single strong emission line, which is identi- 
fied with Lya at z=4.424 ( |Waddington et al. 1999| here- 
after 



W99). VLA J123642+621331 is also identified in 
the ISO supplementary list of sources (Si5^ ~ 23^2/iJy, 
Aussel et al. 1999) but there is no 850 fim sub-mm 
SCUBA detection to a limit of 2 mJy ( [Hughes et al. 1998} ) . 
W99 interpret the source as a dust obscured, nuclear star- 
burst with an embedded AGN. 

The high resolution 0.15" 1.4 GHz VLA-MERLIN ob- 
servations (M99) resolve VLA J123642+621331, show- 
ing that 10% of the flux density resides in an extended 
component lying to the east of an unresolved core. It is 
the core which is detected by the EVN at 1.6 GHz. No 
other sources of radio emission (> 3c) are detected in a 
0.7" x 0.7" field centred on the core. The main question 
is whether the radio emission arises due to star-formation 
processes or AGN activity. The EVN limit on the mea- 
sured core size (< 20 mas) implies a brightness temper- 
ature, T;, > 2 x 10 5 K. While this in itself does not rule 
out a starburst interpretation for the radio emission, the 
high luminosity (~ 10 25 W/Hz i.e. ~ 100 times more lumi- 
nous than the ultra-luminous starburst Arp 220) argues 
strongly for an AGN interpretation for most of the radio 
emission. Further, high resolution, global VLBI observa- 
tions are required to confirm an AGN interpretation for 
this radio source. 

3.3. VLA J123646+621404 

VLA J123646+621404 is a flat spectrum radio 
source located within the HDF itself and is as- 
sociated with a very red, face on, z=0.960 spiral 
galaxy (R98). Rowan- Robinson et al. (1997)| inter- 
pret the ISO 7 /xm emission in terms of a massive 
starburst but the presence of broad Mgll emission 



(Phillips et al. 1997), radio variability (R98) and X-ray 



emission (Hornschcmcicr ct al. 2000) suggest the spiral is 
host to an AGN. We see a 4<r radio source in this field 
that is coincident with the MERLIN position (see Fig. 1). 
We suggest this is compact radio emission from the AGN. 

4. Conclusions 

The results presented here clearly demonstrate the poten- 
tial for VLBI imaging of the micro-Jansky source popula- 
tion. We have reliably detected two radio sources within 
the HDF region. One of these VLA J123644+621133 
shows all the typical properties of a low luminosity, clas- 
sical FR-I radio source. The second, perhaps more inter- 
esting source is VLA J123642+621331. Our detection of a 
relatively compact but high luminosity radio core, further 
supports the suggestion by W99 that this is a dusty, star- 
forming galaxy with an embedded AGN. We also appear 



to have detected radio emission from the AGN associated 
with VLA J123646+621404. 

Global VLBI observations with 3 times the resolu- 
tion and 4 times the sensitivity of the pilot observations 
presented here, are now achievable. With a linear reso- 
lution of ~ 30 pc at cosmological distances, only Global 
VLBI can resolve nuclear star-forming regions like those 
observed in Arp 220 QSmith et al. 1998 ). Targeted obser- 
vations of the optically faint radio source population (of 
which VLA J123642+621331 is just one example) may 
throw new light on the nature of the associated population 
of dusty sub-mm so urces, recently discovered by SCUBA 
( Hughes ct al. 1998 ). In the absence of near-infrared, op- 
tical, ultraviolet and even X-ray emission (due to dust 
obscuration), future VLBI observations will be a crucial 
diagnostic in distinguishing between AGN and starburst 
activity in these systems. 
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